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易 腐 品 多 周期 闭环 物流 网 络 鲁 棒 优 化 模型 
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摘 要 : 针对 不 同 周 期 的 易 腐 品 需求 与 退货 不 确定 性 问题 ， 构 建 了 易 腐 品 多 周期 闭环 物流 网 络 ， 并 设计 了 对 应 的 混 
合 整数 线性 规划 (MILP) 模型 ， 以 实现 最 低 系 统 总 成 本 、 最 佳 设 施 选 址 以 及 最 优 配送 车 辆 运输 路 径 的 决策 。 为 有 效 
规避 不 确定 参数 的 影响 ， 采 用 基 约 束 鲁 棒 方法 ， 将 模型 中 的 部 分 清晰 约束 转换 为 鲁 棒 对 应 式 。 以 上 海 市 果蔬 农产品 
企业 为 实例 ， 通 过 遗传 算法 对 模型 进行 求解 。 结 果 表 明 ， 相 对 单 周 期 而 言 ， 多 周期 系统 具有 动态 性 、 系 统 成 本 更 低 
的 优点 ， 同 时 通过 不 确定 预算 参数 的 变化 分 析 ， 验 证 了 鲁 棒 模型 的 可 行 性 与 有 效 性 ， 进 而 为 不 确定 环境 下 构建 多 周 
期 闭环 物流 网 络 及 降低 系统 成 本 提供 了 借鉴 。 
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Robust optimization model of multi-period closed-loop logistics network for perishable product 


Guo Jianquan, Yang Xiaohua 
(Business School, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: In view of the uncertainty of demand and return for perishable products in different periods, a multi-period model 
of closed-loop logistics network for perishable product and the corresponding mixed integer linear programming (MILP) 
model were proposed, so as to achieve the decision of the minimum total system cost, the best facility location and vehicle 
route. In order to avoid the influence of uncertain parameters effectively, some defined constraints in the model were 
converted to robust corresponding expressions by using cardinality constrained robust approach. Through the case of the 
fruit and vegetable agricultural enterprises of Shanghai, the validity of model was verified by genetic algorithm (GA) . The 
result shows that compared with the single period system, the multi-period system had the advantages of dynamic 
performance and lower system cost, meanwhile, the feasibility and effectiveness of the robust model were verified through 
the analysis of the variation of uncertain budget parameters, which provides a reference for the construction of multi-period 
closed-loop logistics network and the reduction of system cost under uncertain environment. 
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0 3j 并 利用 遗传 算法 求解 。 以 上 文献 尚未 考虑 到 在 参数 不 确定 下 ， 
ES [ 何 同时 统筹 闭环 物流 网 络 设计 与 多 周期 规划 。 
随 着 社会 经 济 的 发 展 以 及 环境 污染 的 加 剧 ， 政 府 和 消费 在 不 确定 环境 下 ， 研 究 闭 环 物流 网 络 ， 更 具有 现实 意义 
(o0 者 对 企业 的 环境 和 社会 责任 提出 了 更 高 的 要 求 避 3。 闭环 物流 。 (4。 大 多 数 研究 采用 随机 规划 、 模 糊 规 划 和 和 鲁 棒 规划 来 开展 
网 络 正 是 顺应 这 一 发 展 趋势 而 产生 的 ， 它 是 科学 高 效 规划 物 ”闭环 物流 网 络 设计 研究 1%1。 然 而 ， 随 机 规划 方法 需要 数据 
流 网 络 的 出 发 点 ， 是 正 向 物流 网 络 与 逆向 物流 网 络 的 集成 ”的 精确 分 布 情况 ， 并 且 目 标 函 数值 会 随 着 场景 数量 增加 而 急 
451。 随 着 闭环 物流 网 络 的 发 展 , 与 其 直接 相关 的 物流 仓储 加 、 — 剧 增 加 ， 由 此 对 模型 计算 提出 了 实质 性 的 挑战 06289， 模 糊 规 
配送 车 辆 路 径 赔 、 选 址 与 分 配 甸 等 研究 的 关注 度 日 益 增加 ， 划 的 不 足 之 处 在 于 ， 其 约束 条 件 难以 权衡 ， 进 而 影响 最 优 决 
并 且 越 来 越 多 的 研究 聚焦 于 系统 性 的 数学 模型 设计 以 及 特定 ，” 策 59。 针对 上 述 规划 的 不 足 ， 不 少 研究 倾向 于 构建 鲁 棒 优 化 
的 框架 分 析 外 。 因 此 本 文 设计 了 易 腐 品 多 周期 半 环 物流 网 络 ”法 09。 重 棒 优 化 能 够 有 效 地 应 对 参数 的 扰动 ， 保 证 物流 网 络 
模型 及 选 址 与 路 径 规划 ， 以 期 为 该 类 问题 提供 系统 性 的 解决 设计 的 稳定 性 。 
思路 ， 并 为 相关 企业 提供 借鉴 。 易 腐 品 是 闭环 物流 网 络 的 重要 研究 对 象 。 文 献 [20] 将 易 
当前 国内 外 学 者 在 该 领域 已 经 开展 了 相关 研究 。 文 献 。” 腐 品 分 为 两 大 类 : 一 类 以 食品 、 蔬 菜 、 人 体 血液 和 使 用 寿命 
[9,10] 通 过 建立 混合 整数 数学 模型 ,来 设计 不 同 产 品 方案 的 闭 。 较 短 的 底片 为 代表 的 易 逝 品 (perishable products); 男 一 类 如 
环 物流 网 络 。 文 献 [11,123] 利 用 元 启发 式 算法 来 优化 多 产品 闭 。 酒精 、 汽 油 以 及 其 他 没有 保质 期 的 放射 性 物质 ， 称 为 衰变 产 
环 供应 链 。 多 周期 性 是 闭环 物流 网 络 的 重要 研究 方向 中。 文 7) (decaying products)。 由 于 易 腐 品 的 保质 期 与 销售 周期 短 、 
献 [1 和 1 研究 了 闭环 物流 网 络 中 具有 序列 相关 设置 的 多 阶段 多 ”退货 品 残 值 低 、 需 求 波动 性 强 等 特性 户 2， 物 流 过 程 的 库存 
产品 多 周期 生产 计划 问题 。 文 献 [15] 为 解决 闭环 物流 网 络 的 。 成 本 、 运输 成 本 很 高 "9。 a Merit 以 更 好 地 
规划 问题 ， 提 出 了 混合 整数 线性 规划 框架 下 的 多 周期 模型 ， 规划 、 控 制 运输 活动 ， 从 而 减少 易 腐 品 在 运输 过 程 的 损失 
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”21。 现 有 易 腐 品 研究 大 多 关注 库存 管理 系统 、 选 址 与 路 径 、 ” 漆 送 至 垃圾 处 理 厂 焚 烧 或 填 埋 1。 
定价 决策 等 方面 es2， 对 于 闭环 物流 网 络 维度 的 探讨 关注 度 易 腐 品 闭环 物流 网 络 网 络 结构 如 图 1 所 示 ， 由 消费 者 自 
不 足 ， 且 缺乏 同时 提 及 不 确定 性 与 多 周期 因素 在 模型 中 的 应 — 提 点 、 易 腐 品 企业 、 中 转 中心 、 易 腐 品 供应 商 、 加 工厂 和 垃 
用 。 圾 处 理 广 六 部 分 组 成 。 

综 上 所 述 ， 针 对 消费 者 需求 和 退货 不 确定 、 易 腐 品 保质 


期 短 的 问题 ， 建 立 了 易 腐 品 多 周期 闭环 物流 网 络 鲁 棒 优化 模 
型 。 为 了 使 企业 实现 最 低 经 济 成 本 目标 ， 本 文 以 最 小 系统 总 


易 腐 品 企业 


成 本 为 目标 函数 ， 包 含 建设 成 本 、 维 持 成 本 、 运 营 成 本 、 信 
息 处 理 成 本 以 及 运输 成 本 。 以 上 海 市 某 果 蔬 农 产品 企业 新 鲜 
直送 水 果 业 务 为 例 ， 通 过 遗传 算法 确定 不 同 周期 不 同 预算 水 和 
平 下 模型 的 适应 度 、 选 址 方案 与 车 辆 路 径 。 以 上 文献 与 本 文 图 1 易 腐 品 闭环 物流 网 络 结构 
的 对 比 情况 如 表 1 所 示 。 Fig.1 Structure of perishable products closed-loop logistics network 
NR m 2 易 腐 品 多 周期 闭环 物流 网 络 设计 
able 1 Literature comparison 
闭环 物流 选 址 与 遗传 不 确定 参 研究 p 2.4 模型 假设 
网 络 KR ”算法 ”数量 WA a) 退货 易 腐 品 只 考虑 第 三 方 再 加 工 与 废物 处 理 两 道 工 
Hiassat A 4 " 易 腐 序 。 
(2017) 品 b) 中 转 中 心 、 供 应 商 、 加 工厂 以 及 垃圾 处 理 广 的 候选 位 
Guo et 4 4 "T 置 与 数量 已 知 ， 并 且 相 同 设施 不 同 节点 的 处 理 能 力 相同 。 
al.(2017) c) 退货 产品 的 加 工厂 与 垃圾 处 理 厂 为 第 三 方 工厂 , 不 考 
Govindan et 4 4 2 打印 ”加权 模糊 数 虑 其 设施 固定 建设 成 本 。 
al.(2017) 机 学 规划 d) 多 周期 情况 下 ， 第 一 阶段 考虑 设施 建设 成 本 ， 其 他 阶 
Zahiri B (2017) ~ 2 医药 ”和 鲁 棒 优 化 段 只 计算 设施 维持 成 本 。 
Guo et F : i "m e) 各 周期 下 设施 的 处 理 能 力 无 差异 。 
al.(2017) A f) 企业 的 信息 处 理 成 本 与 消费 者 的 需求 量 和 退货 量 成 线 
Amin et 4 p 复印 yaspi 性 关系 ， 中 转 中 心 到 供应 商 、 加 工厂 以 及 垃圾 处 理 三 的 运输 
al.(2013) 机 成 本 与 运输 量 和 运输 距离 成 线性 关系 ;中 转 中 心 到 消费 者 
Babaveisi et 是 点 、 自 提 点 SES Iis WA xta ZEE 
J ~ 
al.(2017) 关系 B030。 
Farrokh et 4 i 鲁 棒 模 糊 随 g) 节点 之 间 的 距离 并 非 两 点 间 直 线 距 离 , 而 是 货车 行驶 
al.(2017) 机 规划 的 距离 。 
Hamed et d 22 符号 
al.(2015) i 代表 消费 者 自 提 点 ，iVi ef1,2.…, 了 人 }; j 代表 中 转 中 
m AA vss 4», je(.2...J) 5 s 代表 易 腐 品 供应 商 ，s e{1,2…,5} ; h 代表 
TRGA N 电器 SPAA SEREN, hel2H}; p 代表 加 工厂 ，p ef12..,P} ;1 
李 伯 棠 等 4 i 发 动 ibid 代表 垃圾 处 理 厂 ，1e{12.…,L} ; t RRA, £602...) ; v 
(2017) 机 代表 运输 车 辆 ，v eV ; R 代表 运输 路 线 ，R =12,…,R o 
、 易 腐 。 23 参数 
B E X Ee us. t0 Fr 、 Fi 分 别 为 第 1 周期 ， 各 节点 设施 建设 的 固定 成 本 ; 
1 ”问题 描述 maj, 、ma; 分 别 为 第 1 周期 , 各 节点 设施 的 维持 成 本 ; opt. ovi 
EL 分 别 为 第 t 周期， 各 节点 的 运营 成 本 ; K 代表 单位 车 辆 最 大 
单 周 期 物流 配送 中 ， 忽 略 了 消费 者 对 易 腐 品 的 需求 与 退 ” 运载 量 ; wi 代表 易 腐 品 企业 的 单位 信息 处 理 成 本 ; is 、 


货 是 动态 变化 的 因素 ， 从 而 导致 闭环 物流 网 络 系统 中 设施 长 。 wn U ts utp wt 分别 代 表 两 点 之 间 的 单位 运输 成 本 ， 
期 选 址 在 实际 决策 中 出 现 偏差 。 而 多 周期 物流 配送 ， 充 分 体 diss dy dis. dines dis, 分 别 代 表 两 点 之 间 的 距离 ， ca ca) 、 
时 配送 路 径 的 灵活 性 ， 并 且 有 效 解决 不 同 周 期 参数 值 变化 对 co、 ca 分 别 代表 各 节点 的 最 大 处 理 能 力 ;， an 代表 回收 产品 的 
设施 长 期 选 址 与 数量 的 影响 ， 使 得 研究 更 贴切 实际 。 在 多 周 MIMIX: i 代表 消费 者 自 提 点 第 ANERE TREN 
关系 统 中 ， 根 据 产品 易 腐 程度 、 市 场 环境 以 及 易 腐 品 企业 运 。 “ 费 者 自 提 点 第 1 周期 退货 量 。 
营 的 实际 变化 等 因素 进行 周期 划分 ， 并 利用 计算 得 出 的 最 优 24 决策 变量 

设施 选 址 安排 进行 各 个 周期 的 物流 配送 作业 。 站 为 O1 变量 ， 若 在 第 周期 选择 中 转 中 心 ， 则 
在 某 一 周期 内 ， 中 转 中 心 配送 仓 作为 闭环 物流 网 络 的 分 as snl UEM REL 

销 商 ， 需 对 从 供应 商 采 购 的 易 腐 品 进行 分 类 、 包 装 、 贴 签 等 和 
作业 。 配 送 运输 过 程 中 ， 乘 坐车 辆 y 的 工作 人 员 将 产品 按 企 Li ano OIAR, EAA. IERRA 
业 预 售 的 订单 需求 配送 至 正 向 物流 终端 (消费 者 自 提 点 ), 与 7 0 UO 8 0-1 变量 ， 若 在 第 + 周期 选择 候选 垃圾 
此 同时 ， 完 成 提取 退货 品 的 一 系列 活动 。 在 逆向 物流 中 ， 中 A, Y= AUY =0 下 为 0-1 变量 ， 若 车 辆 v 在 t 周 
转 中 心 回收 仓 对 各 自 提 点 的 退货 易 腐 品 分 类 处 理 ， 将 DES R 条 路 线 运输 易 腐 品 ， 则 只 =1 否则 Y# -0; YE, YE, YA 
定 残 值 的 产品 送 至 加 工厂 进行 再 循环 工序 ， 后 续 把 易 腐 品 残 。 为 0-1 变量 ， 若 在 t 周期 第 R 路 线 上 ， 分 别 将 产品 从 中 转 中 
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MES, F: 易 腐 品 多 周期 闻 环 物流 网 络 鲁 棒 优化 模型 


心 运输 到 


Ea HIER i 运输 到 自 提 点 ?， 以 及 从 自 提 点 


运输 到 中 车 中 心 ， 则 Yi =l ` YSL, Yi =l, 否则 Yi -0. Yj-0 ` 


vg =0 


X 


3 


t 
Tux 


Xy. Xp. Xy. RIRA t JALULIR A [e] x dp E e 


2.5 MILP 模型 的 建立 


Wsin = 


单 周 期 的 固定 建设 成 本 与 设施 维持 成 本 为 
YR 5377 ^Y Y Yx ma! HY XY masses, +x) (1) 


多 局 


Ww 
ZH 
zn 


期 的 固定 建设 成 本 为 


I J 
Wmi (t=1)= X F + YF! (2.1) 
i=l j=l 


期 的 设施 维持 成 本 为 


mE Wd 
w, (22) - Yr YI mai D Y Yima; * (x +x) (2.2) 
j=l i=l jaia 


司 期 的 运营 成 本 、 信 息 处 理 成 本 及 运输 成 本 为 
m= YYYYons, 十 » ; yyron (x 9x5) 
F n E £ d » J S 
SSe eee Eines 
AL jal i0 tal j=l s20 


dis jut; + yyy yy Yixidis put ip + )»5»»37 YY; x dis ut, + 


t=1 j=1 p=0 i=l 1=1 j=l p=0 1-0 


Y d Gg T Ir rA P: 
SSY IHY, Y Yt idis, 3-3 9 YYP dis, 
veV t=l i=0 i'=0veV t=l i=0 j-0veV 


t=1 j=l i=0 
(3) 
周期 系统 总 成 本 为 ”min Zs =W +W (4) 
多 周期 系统 总 成 本 为 ”minZ = wo +m 6) 
s. t. 
T Pr 
YYXUORYozi,vi: (6) 
t=l Pl 
Ye=0, Vit (7) 
f rf 
> YY =i, vut (8) 
i-0 7=0 
x; ZdiY$ , Vi j,R,t (9) 
CREE i 
»»ES = Xy | Vs f (10) 
mu 
Xij 2 Yg V? j,R,t (11) 
MS i J J 
a=’, /a= (a). Yp b t (12) 
i=l j=l jal j=l 
YR <l, YR, t 
> (13) 
I J* 
TYWYMEREK,VRbovet (14) 
i-0 i'-0 
xytxpXca, Vi jt (15) 
I J 
a? V uj Sca, Vp. t (16) 
ja 
tog 
0-009 Xj €cap Vlt (17) 
il j=l 
r EI 1 J 
22295 2.2/4 sca, vt (18) 
j=l i i=l j=l 
J 
rl vt (19) 
jal 
b 
XY: 21, vt (20) 
pz 
L 
XYL vi (21) 
l=1 
yr>1, Vt 
22 Q2) 
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Y, Y» Yb YR, Y Yi,Y$e(0lp Vi b p L v R Q3) 
Mp Xj Xp Xj x,20 Integer, Vs. j i p l (24) 
模型 中 ， 式 (60 (7) 表示 两 个 自 提 点 间 仅 有 一 条 线路 ; 


X (D 表示 一 个 自 提 点 只 属于 一 条 路 线 ; SX (9) ~ (12) 
表示 闭环 物流 网 络 中 ， 易 腐 品 产品 流量 平衡 关系 ; X CIS) 
(14). 表示 每 条 线路 仅 有 一 辆 车 配送 ， 并 且 每 条 线路 运输 总 
量 不 超过 车 辆 运载 量 ; 式 (15) ~ (18) 分 别 表 示 各 个 设施 
的 容量 约束 ; 式 (19) ~ 21) 表示 至 少 选 择 一 个 中 转 中 心 、 
一 个 食品 加 工厂 和 一 个 废物 处 理 厂 ， 式 (22) 表示 至 少 有 
条 线路 完成 整个 闭环 供应 链 系统 的 配送 ， 式 〈23 ) 表示 变量 
为 0-1 变量 ， 式 (24) 表示 运输 量 为 非 负 整数 。 


3 ”多 周期 下 鲁 棒 优化 模型 设计 


在 闭环 物流 网 络 问题 中 ， 考 虑 易 腐 品 不 同 周 期 、 需 求 与 
退货 不 确定 因素 对 选 址 、 路 径 决 策 至 关 重 要 。 为 有 效 规 避 不 
确定 参数 的 影响 ， 采 用 基 约 束 鲁 棒 方 法 ， 将 模型 中 的 部 分 清 
晰 约束 转换 为 鲁 棒 对 应 式 ， 并 结合 对 偶 理论 ， 对 非 凸 非 线性 
的 不 确定 鲁 棒 规划 进行 线性 变换 0%19。 

假设 每 个 节点 的 需求 量 4 取 值 于 盒子 (box) 不 确定 集合 ， 


即 有 界 的 对 称 区 间 集 合 。 必 <[d -ww arm]. KP d deg 


t 周期 需求 量 的 常规 值 ， 的 最 大 偏离 


值 。 在 约束 式 (9〉 中 ， 引 入 不 确定 预算 参数 ， 来 调节 结 
果 的 保守 性 ， 其 中 Tie [0, m], m 表示 每 个 周期 下 需求 点 
不 确定 的 个 数 。 在 现实 的 易 腐 品 闭环 供应 链 中 ， 并 非 全 部 的 
需求 量 都 是 不 确定 值 。 本 文 构建 的 鲁 棒 优化 方法 是 为 了 保证 


最 多 有 蕊 个 4d; 在 其 区 间 内 波动 ， 对 应 的 需求 量 扰 动 为 


d 表示 第 t 周期 需求 量 


(ri-rs)d 。 由 此 ， 式 (9) 转换 为 
x, di Yk +z(di Ti), Y ij, Rt; 
* di- tYR d [wd tYR 

其 中 : z(d; S gent amp Hu T; JLYR)} (25) 
将 z (4 T1) 进行 线性 转换 : 
I 

z(d,',V$) - max(diY4ó/ : > 6 STY,0< 6 <1, Viel (26) 
i=l 

其 中 : * 表示 对 数 * 向 下 取 整 ，S eR EISE dE z(dy. T) 


为 保护 函数 ， 6? 为 辅助 变量 。 

当 I=0， 此 时 表示 所 有 自 提 点 的 需求 量 都 均 未 偏离 均 
值 ， 这 种 情况 下 ， 重 棒 问 题 变 为 确定 问题 ; T =m 此 时 所 
有 需求 值 都 偏离 均值 , 取 值 于 区 间 最 大 偏离 值 , 在 此 情况 下 ， 
所 得 的 物流 成 本 过 于 保守 。 企业 可 以 通过 TT; 值 来 调节 目标 值 
的 鲁 棒 性 。 由 于 式 (C250 有 界 ， 根 据 对 偶 定 理 ， 其 对 偶 问 题 
X QD 有 界 且 所 得 结果 一 致 。 


I 
uri - Y vtiutevi 2dYR, uf>0, vi20, v1, j; R, |] (27) 
i=l 


min 


根据 式 (25) ~ (27)， 需 求 量 约束 式 (9) 的 鲁 棒 转换 
等 价 式 如 下 : 


I 
x, 2 dYE +u ri +Ý Vv, vij. R t (28) 
i=l 
u; +v; ZdiYS , Vij, Ri (29) 
u20 20, vi jt (30) 


同 理 ， 将 消费 者 的 退货 量 约束 式 〈11) 进行 鲁 棒 对 应 式 
优化 转换 : 
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V(R) 


H FA 


su 


PPERBPEBRDPEBRBPET] 


" (R) Dg 


1 1 2 48 


(b) 染色 体 组 2 
(b) Chromosome 2 


SET 


(c) 染色 


色 体 组 3 


(c) Chromosome 3 
图 2 ”染色体 组 


Fig. 2 Chromosomes 


R2 遗传 算法 参数 


Table 2 Parameter of GA 
种 群 数量 交叉 率 变异 率 和 迭代 周期 
数值 20 0.8 0.25 100 
b) 适应 度 评 估 与 选择 。 染色体 的 适应 度 反 映 了 被 选择 概 


率 。 
利用 


收敛 性 


NES 


更 高 的 适应 1 


在 满足 式 (6) ~ (8) (10) (12~24) (28~33) 情况 下 ， 
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x eU eerie vj. R t G1) pend Men ONU esci edid 
d) 终止 条 件 。 若 达到 预先 设 定 的 最 大 进化 代数 ,终止 算 
esq ym, Vh t G2) 法 ， 否 则 转向 步骤 b)。 
e, 20, q20,vVijt (33) 5 算 例 
5.1 数据 来 源 
4 中 传 算法 编写 过 程 以 上 海 市 某 果 蔬 农 产品 企业 新 鲜 直送 水 果 业 务 为 例 ， 假 
遗传 算法 CGA) 是 基于 遗传 学 机 理 和 达尔 文 自然 选择 理 设 配 送 车 辆 运输 水 果 产 品 单 周 期 长 度 为 一 周 ， 多 周期 中 一 周 
论 的 元 启发 式 优化 技术 中 ， 适 用 于 搜索 近似 最 优 解 24。 分 为 三 个 周期 (精确 的 周期 长 度 划 分 以 产品 易 腐 程度 、 市 场 
有 内 在 的 隐 并 行 性 、 全 局 寻 优 能 力 、 自 动 调整 搜索 方向 等 特 环境 以 及 企业 运营 的 实际 变化 为 准 )。 通 过 对 企业 调研 ， 三 个 
点 ， 被 广泛 用 于 解决 大 规模 运输 规划 、 物 流 网 络 构建 、 配 送 ”候选 的 中 转 中 心 坐标 分 别 为 (5.4，5.0)、(5.6，3.8)、(7.6， 
车 辆 路 径 等 问题 B133。 在 此 基础 上 ， 本 文采 用 GA 来 解决 进 — 370: 两 个 候选 的 加 工厂 坐标 为 (1.6，3.9)、(1.2，5.1); 2 
行 多 周期 闭环 物流 网 络 鲁 棒 设 计 、 选 址 及 路 径 规划 问题 。 本 个 候选 的 垃圾 处 理 厂 坐 标 为 (2.8，7.1)、(1.5，7.9); 一 个 
文具 体操 作 如 下 所 示 : 选择 的 水 果 供 应 商 坐 标 为 (2.8，4.4); 八 个 选择 的 消费 者 水 
a) 染色 体 编 码 与 初始 化 。 本 文 的 染色 体 (J+P+L+ V(R) 果 购 买 /退货 点 坐标 及 各 周期 的 需求 与 退货 预测 值 如 表 3 所 
+JI(R)+II(R)+IJ(R)+SJ+JI+HIT+IJ+JP+POŚX H. K 2(a)(b) 示 ， 其 他 参数 数值 见 表 4。 
的 染色 体 为 二 进 制 串 (1 或 0), 表示 设施 J P. L 是 否 建立 、 表 3 消费 者 自 提 点 坐标 、 各 周期 需求 量 及 退货 量 
车 辆 V 是 否 在 路 线 R 上 运输 以 及 路 线 R 上 节点 是 否 有 被 选 Table 3 Pick-up point coordinates, demanded and returned quantity 
Tí; 图 2(c) 的 染色 体 对 应 的 数值 表示 实例 中 节点 间 的 易 腐 品 received at pick-up point in each period 
运输 数量 。 在 GA 的 过 程 中 ， 自 动 确定 消费 者 自 提 点 、 中 转 (X, Y) d db db mg m m 
中 心 、 加 工厂 及 垃圾 处 理 厂 是 否 开放 以 及 规划 相应 路 线 。 本 L (13.6, 10.1) 53 57 50 09 09 0S8 
文 的 GA 参数 如 表 2 所 示 。 L, (2.7, 12.0) 56 53 59 08 07 0.8 
A (10.1, 3.2) 74 78 73 12 12 14 
L (1.5, 2.2) 60 62 62 10 10 10 
FREESESEEEEU nh 67,69 58 50 52 10 08 09 
(a) 染色 体 组 1 I, (10.4, 7.3) 69 65 65 12 1l 1.1 
(a) Chromosome 1 L (10.2, 1.4) 7.8 73 73. 12. 11 13 
I (8.1, 4.9) 57 57 58 10 09 09 


ARA 其 他 参数 相关 数据 


Table 4 Relevant parameter data 


参数 取 值 参数 取 值 ”参数 取 值 

F; 6000 元 ca, 10t Vo 3 辆 

Fj 30000 元 ca, ]00t ui, 20 元 -t7 

ma 1.5 元 :ta ca, 18t utj 60 元 -km-! 

ma' 4 jb t^! ca,  ]Ot ut; 376: (kmt ) 

opi 10 元 t^! K  100t utj -ut;; 80 元 km! 

op; 50 元 t^ a 70% ut; —ut, 4 元 (km-t ) 

为 考虑 水 果 需 求 量 与 退货 量 不 同 鲁 棒 水 平 下 对 最 优 策略 


的 影响 ， 分 别 对 Dir 赋值 0、3、6。 在 此 基础 上 ， 将 两 类 
不 确定 预算 交叉 产生 的 九 种 情况 


casel (F4 20, T; =0) 为 确定 情况 , case9 (T4 26. T; =6) 为 最 


情况 。 假 设 需求 与 退货 的 最 大 偏 


带 入 模型 中 ， 其 中 
保守 
5 值 为 常规 值 的 15%， 即 


d =15%xd 、rm =15%xr ， 那 么 


离 值 也 越 大 。 
5.2 运算 结果 


通过 软件 


算出 的 目标 函数 值 结果 与 设施 选 址 方案 如 
a) 观测 九 种 
随 着 不 确定 预算 TL 与 ;的 增加 ， 每 个 周 


从 鲁 棒 性 角度 分 析 可 得 : 


选择 操作 避免 遗漏 ， 


目标 函数 (1) ~ C50 计算 各 个 个 体 的 适应 度 值 ， 适 
度 值 可 直观 反映 个 体 的 优 劣 。 


< 用 轮 盘 赌 策 


出 适应 度 大 的 个 体 作为 父 代 ， 提 
c) 交叉 与 变异 。 


生 。 变 异 使 


交叉 使 子 代 同时 继承 父 


攻略 进行 选择 。 通 过 评估 适应 度 值 ， 
高 全 局 收敛 性 。 
母 的 基因 


HEr 
JE In] E 


某 个 基因 以 一 定 概率 随机 发 生 突变 。 


好 的 鲁 棒 怕 


自 提 点 需求 /退货 量 越 大 ， 偏 


过 MATLAB 2016a 对 遗传 算法 


进行 编码 , 迭代 运 
K 5 Wn. 
情况 的 纵向 比较 ， 


期 的 需求 与 退货 扰 


动量 均 呈 上 升 趋势 ， 系 统 总 成 本 随 之 递增 。 上 述 结 果 表 明 ， 
当 考 虑 的 不 确定 量 越 多 时 ， 为 确保 闭环 物流 网 络 系统 具 
生 ， 整 个 系统 会 增加 设施 点 的 开 


和 具有 更 


放 ， 同 时 系统 总 成 


本 上 升 ; b) K 


3 中 ， 在 不 考虑 周期 因 
的 各 个 阶段 变动 趋势 相近 ), 退货 量 


EN 


〈 单 周期 与 多 周期 


仅 为 需求 量 的 1596-2096, 
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忆 此 相同 的 不 确定 预算 变化 下 ， 需 求 扰动 量变 化 更 大 ， 使 得 腐 品 闭环 物流 网 络 ,表明 了 和 鲁 棒 优化 模型 的 可 行 性 与 有 效 性 。 


其 鲁 棒 性 影响 更 大 ， 成 本 的 上 升 趋势 更 显著 ， 如 case3-4、 综 上 ， 和 鲁 棒 优化 模型 考虑 了 参数 不 确定 性 下 的 闭环 物流 网 络 
case6-7 变动 情况 所 示 ; c) 即使 在 需求 量 、 退 货 量 和 周期 都 。 的 运作 ， 同 时 模型 鲁 棒 性 与 系统 总 成 本 之 间 存 在 权衡 关系 : 


发 生变 化 的 时 候 ， 重 棒 优 化 后 的 解 依然 可 以 有 效 地 应 用 于 易 ” 模型 鲁 棒 性 越 强 ， 系 统 总 成 本 越 大 。 
表 5 不 同 鲁 棒 水 平 的 各 周期 GA 适应 度 及 选 址 情况 


Table 5 GA fitness and location in each period under different robust levels 


多 周期 "m 
第 一 周期 第 二 周期 第 三 周期 多 周期 gap 
成 本 选 址 成 本 选 址 成 本 选 址 成 本 和 成 本 选 址 
casel(T7 =0, I7 =0) 879038 — jap,lh — 10063  Jjopslh — 8295 — j.p. 106296.00 1309038 joj» p.h 0.23 
case2(T¥ 20, Iz =3) 88486 — j.p. 10613  j.p:h 8846 ,pl 107944.81 131625 jj» Puoh 0.22 
case3(T4 =0, Tv —6) 89335  jL.p.l 11463  Jjapsh 9696 — pi.pspsjpP.h 110493.53 132613 0.20 
case4(T4 23, Iv =0) 91420 — j.p 13567  Jj.peh — 11800 — j.p 116787.05 135312 jojo Pi Pəh 0.16 
case5(T4 =3, I7 =3) 91969 j Poh 14117 j.peh — 12350 — j.p 118435.77 135999 六 万 ,PP 0.15 
case6(TY 23, I7 =6) 92317 jo Poh 14467  J.Psh — 12700 — j.p. 119484.49 136486 | Jojo Pi» D» 0.14 
case7(V4 26, T; 20) 934030  j.p.h — 15571 — jepeh 13804 — jsp 122778.0 | 138186 — Jj. jz; Pi Pəh 0.13 
case8(T4 =6, I =3) 93751 jpuh — 15922 jpuh — 14154 — jo pil 123826.73 138673 jù Jz Pep 0.12 
case9(T4 26, T; =6) 94199  j.p.l 16372 j.p] — 14604 — j.p. 125175.45 139260 — Jj. js Po D. 0.11 
3E: gap=(| 单 周期 成 本 -多 周期 成 本 和 |)/ 多 周期 成 本 和 。 
160000 
140000 + 
120000 
4m 100000 一 一 单 周 其 
E 第 一 周期 
Ej 80000 
e —— 第 二 周期 
[i 60000 >- 第 三 周期 
40000 
20000 
和 
0 
1 2 3 4 5 6 7 8 9 
场景 
图 3 不 同 鲁 棒 水 平 的 总 成 本 变化 情况 j RS i 
Fig.3 Total cost variation at different robust levels ala arah pm LACER . D TIONS 
从 多 周期 角度 分 析 可 得 : a) 横向 观测 各 个 情况 不 同 周期 图 5 第 一 周期 路 径 图 (case 9) 
不 同 阶 段 的 系统 总 成 本 与 选 址 情况 ，GA 和 迭代 结果 表明 ， 多 Fig.5 The first period path diagram (case 9) 
周期 总 成 本 更 低 ， 设 施 选 址 动态 变化 ， 更 加 符合 实际 情况 ， E y ) R ET 
但 由 于 考虑 的 周期 长 度 较 小 ， 设 施 的 固定 建设 成 本 占 系 统 总 ? 


成 本 比重 大 ， 所 以 单 周期 与 多 周期 总 和 差 值 较 大 ， 差 值 比 在 
0.11~0.23 之 间 ; b) 随 着 不 确定 预算 TC. 的 增 大 ， 系 统 总 
成 本 增加 ， 固 定 建设 成 本 比重 减少 ， 单 周期 与 多 周期 总 和 但 
差 值 递 减 ; c) 此 外 ， 如 图 4-7 所 示 ， 在 最 保守 的 情况 下 
(case9)， 单 周期 的 设施 数量 以 及 配送 车 辆 高 于 多 周期 各 个 
阶段 ， 同 时 ， 多 周期 的 三 个 阶段 设施 节点 选 址 虽然 一 致 ， 但 
配送 车 辆 运输 路 线 根 据 消 费 者 实际 需求 变动 而 作出 相应 调 
节 ， 更 加 符合 实际 情况 ， 由 此 说 明 多 周期 下 闭环 物流 网 络 的 


设施 选 址 及 路 径 更 加 具有 灵活 性 ， 并 且 系 统 总 成 本 更 低 。 (case 9) 
K er Y xc. N Fig. 6 The second period path diagram (case 9) 
m. K xm—e U€U50 vag Y 
x w s 
-À Xe 


图 4 单 周 期 路 径 图 (case 9) 图 7 第 三 周期 路 径 图 (case 9) 
Fig.4 Single period path diagram (case 9) Fig.7 The third period path diagram (case 9) 


录用 定稿 mta, # 
6 ”结束 语 


本 文 将 易 腐 品 设 为 研究 对 象 ， 以 最 低 系 统 总 成 本 、 最 佳 
设施 选 址 以 及 最 优 运 输 路 径 为 出 发 点 ， 建 立 了 多 周期 闭环 物 
流 网 络 鲁 棒 优化 模型 ， 并 利用 遗传 算法 求解 模型 。 结 果 分 析 
表明 ， 重 棒 优 化 模型 能 够 有 效 地 控制 易 腐 品 闭 环 物流 网 络 的 
需求 与 退货 不 确定 性 ， 同 时 多 周期 系统 的 设计 更 具 灵 活性 、 
降 成 本 的 优点 。 在 此 基础 上 ， 增 加 回收 工序 、 假 设 、 约 束 等 
条 件 ， 定 制 出 的 模型 也 可 适用 于 玻璃 容器 、 包 装 纸 盒 、 服 包 
等 其 他 可 循环 产品 的 多 周期 闭环 物流 网 络 鲁 棒 设计 及 路 径 优 
化 。 


zi 


在 以 往 研究 基础 上 考虑 了 易 腐 品 闭环 物流 网 络 的 多 周期 
及 不 确定 因素 ， 使 得 本 文 模型 更 具有 现实 指导 意义 。 然 而 ， 
本 文 在 模型 设计 时 提出 了 多 条 假设 条 件 ， 在 后 续 研 究 中 ， 可 
适当 放松 部 分 假设 使 得 模型 更 加 贴近 企业 实际 ， 另 外 ， 在 本 
文 的 计算 案例 里 ， 仅 考虑 了 时 长 为 一 周 的 多 周期 〈 多 次 ) 配 


E 


送 问 题 ， 甚 中 周期 期 限 与 次 数 选取 代表 性 尚 有 不 足 ， 下 一 步 
研究 有 必要 扩充 其 适用 性 与 代表 性 ; 最后， 解决 闭环 物流 网 


络 不 确定 问题 的 方法 不 止 本 文 提 及 的 鲁 棒 优 化 法 ， 如 何 将 其 
他 不 确定 方法 〔 如 模糊 约束 法 、 随 机 规划 法 、 鲁 棒 优 化 法 等 ) 
应 用 到 模型 中 并 分 析 各 自 优 劣 势 ， 还 待 进一步 探讨 。 
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